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ABSTRACT 

Finite-amplitude    standing   waves    in   air    at    ambient    temperature 
contained    in   a   rigid-walled  rectangular    cavity   were    experimentally 
investigated.      The  pressure   waveform    in    the   cavity   was    analyzed   for 
harmonic    content   and   compared   to    the    theory    of   Coppens    and   Sanders 
as   modified   to    include    empirical    losses.      The   waveforms    for    the 
nOO   and  OnO    configurations    were   analyzed  at    strength   parameters    of 
about   0.5,    0.75,    and    1.0,    where    the   strength   parameter    is    defined  as 
MQ    (y    +    l)/2.       M    is    the   Mach    number    of    the    fundamental,    Q      is    the 
quality    factor    of   the    fundamental    resonance,    and  y    is    the    ratio   of 
specific   heats.       It    is    found   that    the   theoretical   model   accurately 
predicts    the   shape   and  approximates    the   amplitudes    of    the   harmonic 
components    of   the   finite-amplitude   standing   waves. 
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I.   INTRODUCTION 


The  purpose  of  this  research  is  to  make  a  detailed  investigation 
of  finite-amplitude  standing  waves  in  air  at  ambient  temperatures  in 
a  rigid  walled  cavity,  for  comparison  with  the  theory  of  Coppens  and 
Sanders  [l]  modified  to  include  observed  energy  losses.  This  in- 
vestigation is  a  continuation  of  the  experimental  works  of  Beech  [2] 
and  Winn  L3j,  with  an  extension  to  a  rectangular  cavity  vice  a  long 
round  tube. 

A  rectangular  cavity  was  chosen  as  it  was  a  step  toward  the 
study  of  two  and  three  dimensional  waves,  although  this  investi- 
gation concerns  only  one  dimensional  standing  waves.   In  addition, 
it  was  thought  that  a  rectangular  geometry  would  produce  over- 
tones which  would  be  less  harmonic  than  those  for  the  more  sym- 
metric geometry  of  the  round  tube.   The  dimensions  of  the  cavity 
were  chosen  in  an  effort  to  space  out  the  system  resonances  to 
avoid  degeneracy  or  near -degeneracy .   This  was  accomplished  with 
some  success,  and  thus  allowed  for  the  accurate  location  and  anal- 
ysis of  the  first  few  overtones  of  each  of  the  two  modes  studied. 


II.       BACKGROUND   AND  THEORY 

A   finite  amplitude   wave    is   a    solution    of   the   non-linear    wave 
equation    that    describes    the  real-life   environment    that    confronts 
large   amplitude  acoustic    waves.      This    finite-amplitude   wave   re- 
quires   second   order    terms,    or    higher,    to    describe    the   waveform. 

It    has    been    known    for    many    years    that    hydrodynamic    equations 
describing    traveling   waves    in   a    non-diss ipat ive   medium  predict    a 
change    in    the   waveform.      Kirchoff    L4J    and  Lamb  [5]    were   among   the 
early    investigators    in    this    area.       In   more   recent    years    there   have 
been   many    theories    and   methods    of    solution   proposed   to    describe 
this    distorted   traveling   waveform    in   a    non-dissipative   medium. 
Amoung   these   are    the   works    of   Fay    [6],    and    Blackstock   [7]    who  pro- 
duced an    extensive   theoretical    treatment    of    low  amplitude   non- 
linear   waves,    and  made    use   of    Fourier    and  power    series 
approximations . 

When   the  problem  was    extended   to   diss ipat ive   mediums,    it   was 
shown    very    early    [8]    that    waves    in   absorptive   media    exhibit    finite- 
amplitude   characteristics    only   at    high    strengths,    different    from 
waves    in   a   non-dissipative   media.      More   recently,    an    important    con- 
tribution   to   the   understanding    of   finite   amplitude   waves    in   a    dis- 
sipative  medium   was   presented   by    Fay    LQ3.      He,    similar    to   a    later 
work   by    Blackstock   [lOj,    formulated   a    Fourier    series    solution    to 
describe   the    finite-amplitude   waveform    for    a   plane   wave    of    infinite 
extent.      More   recent    theoretical   approaches    to    finite-amplitude 
effects    have    been    of   two   basic    types:    an   approximation   solution, 


such  as  that  proposed  by  Fox  and  Wallace  Cll],  and  a  perturbation 
solution  like  that  of  Fay  C9l.   The  approximation  solution  is  based 
on  assumed  attenuation  mechanisms  and  manner  of  overtone  population. 
The  perturbation  solution  makes  use  of  perturbation  methods  and 
Fourier  series  representation  of  the  waveform  (.10,12,13],   A  per- 
turbation analysis  which  was  extended  to  the  sixth  order  by  Keck 
and  Beyer  C 1 3 ]  is  of  particular  importance  to  this  investigation, 
as  will  be  evident  later  .   The  solution  obtained  by  Keck  and  Beyer 
is  valid  near  the  source,  unlike  previous  perturbation  solutions 
of  finite-amplitude  traveling  waves  in  dissipative  media. 

Although  most  investigations  have  been  confined  to  traveling 
waves,  there  has  been  some  investigations  of  standing  waves.   For 
example,  Keller  L l4l  extracted  solutions  for  standing  waves  in  a 
closed  tube.   In  that  investigation,  Keller  assumed  a  non-dissi- 
pative  medium  which  resulted  in  infinite  amplitude  waves  at  re- 
sonance, yet  still  obtained  useful  results  for  frequencies  close 
to  those  of  resonance.   Coppens  and  Sanders  LIJ  developed  an  ex- 
tension of  the  Keek-Beyer  perturbation  approach,  that  include  wall 
losses  as  predicted  by  the  Rayleigh-Kirchoff  theory  [l5l.   This 
approach  was  carried  to  the  tenth  order  by  Ruff  [l6],  at  which 
point  it  was  found  that  the  mathematics  become  divergent. 
Coppens  [l7l  later  expanded  the  theory  and  used  a  Fourier  decom- 
position technique,  which  both  Beech  L  3*1  and  Winn  L  2  J  experimentally 
investigated.   It  was  found  that  the  tubes  that  they  investigated 
did  not  behave  in  the  manner  predicted  by  Rayleigh-Kirchoff  loss 
mechanisms.   This  investigation  follows  up  the  observation  that 


the    tubes    exhibited   non-Rayleigh-Kir chof f   behavior,    by    testing    the 
theory   using   empirically   determined   losses. 

In    the   theory    of   Coppens   and  Sanders,    a    one    dimensional,    non- 
linear   wave    equation   with    dissipative   terras    which    correspond    to 
those   encountered      in    the   discription    of   plane   standing   waves    in   a 
rigid-walled  cavity,    is    developed.      The   wave   equation    in   this 
theory    can   be   Fourier    decomposed  to   take    the   form 

l&  -*$•■}%■.  m° 

v  o 
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and   contains   a    dispersive    term   and  a    dissipative    term,    both   based 
on    the   Rayleigh-Kir choff   theory    of  wall    losses    [l5l>    and 
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CO   and  CO       are  the  half  power  frequencies  above  and  below  the 
fundamental . 

In  the  description  of  waves  in  cavities,  Eq.  1  can  be  refor- 
mulated Ll8j  to  introduce  the  actual  cavity  resonances  and 
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generalized  to  include  an  arbitrary,  empirical  absorption  coeffi- 
cient for  each  mode.   The  result  is 
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acoustic  pressure 
density  of  medium 
thermodynamic    speed   of   sound   in   medium 


=    mach   number 
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P..    =   peak  pressure    of   fundamental    component    of   wave 

R      =    Fourier    coefficient    of   n         harmonic    component,    normalized 
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0      =   phase  angle    of   n        harmonic   component,    where   0=0 


This  equation  can  be  manipulated  into  a  general,  time-independent 
form  which,  when  evaluated  at  the  rigid  boundaries  of  the  cavity, 
can  be  further  simplified  to  produce  the  following  form  amenable  to 

computer  calculation: 

n  =  l 

H  R  COS  (0    -0    )     =   MbQ,    ^[^S  R.R        tftf(0.+0        .) 

n   n       v    n      n'  1        i        j=x      J    n-j  j  n-J 

CO 

-  E  R        .R  COS  ( 0        .-0.) 

j  =  1      n+j    j      K    n+j        V 

and   a    similar    form   with    "sin"    in   place    of    "  COS  " , 
where 
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and 
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tan    G       =  "    1       "r 
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As    can   be   seen    from   the   above    equations,    it    is    important    to   know 

the   Q(N) 's    and  CO      's    of    the   system    for    the    input    of   the    theoretical 

'         nr  '  r 

predictions  calculated  by  the  computer .   It  is  convenient  to  de- 
fine a  quantity,  E(n),  that  indicates  the  position  of  CO   relative 
to  the  classical  harmonic  frequencies,  nCO   .   E(n)  is  defined  by 

E(n)  =  (CO  /(nco.  ))  -  1. 
\  /     v  nr/  \    lr  /  / 

The  values  of  Q(n)  and  E(n)  are  directly  determined  from  the 
infinitesimal  amplitude  analysis  of  the  cavity  and  are  the  input 
quantities  for  the  computer  calculations. 
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III.   EXPERIMENTAL  CONSIDERATIONS 

A.   APPARATUS 

A   block   diagram   of   the   experimental   apparatus    is    shown    in 
Figure    1.      Finite-amplitude   waves    were   generated    in   air    at    ambient 
temperatures    in   a    rectangular    cavity   with    interior    dimensions    of 
6.97    cm.    by   20.96    cm.    by    3O.28   cm.      The   walls    of    the   cavity   were 
constructed   of   O.985    in.    thick  aluminum   plates,    milled  and    screwed 
together    to  provide   tight,    right   angle    joints.      These   joints    sealed 
by   applying    thin    layers    of   silicon   grease   to  adjoining   surfaces 
before    screwing    them   together.      Acoustic    waves    were   generated    in 
the    cavity   by    exciting    the   air    with   a   plane-faced  piston.      The 
piston    face   gained  access    to    the    interior    of   the   cavity   by    way    of 
an    orifice  bored    in    one   of    the  cavity   walls.      The   piston   was    fitted 
with    two  0-rings    and   the  piston   orifice   was   bored   such   that    the 
piston   with   0-rings    would   make   an   airtight    fit.       It   was    found   that 
the  piston   motion   was    restricted  by    the    second  0-ring;    therefore 
the  piston   was    inserted   into    the   orifice    only   as    far    as    the    first 
0-ring.      Lubrication   and   further    sealing    of   the   piston-orifice 
interface   was   provided   by    the   same   high    vacuum   silicon    grease   used 
to   seal    the   wall    joints.       The   piston   was    attached    to  an   M-B   Elec- 
tronics  Model    EA    1500   Exciter.      This    exciter    was    powered  by    two 
M-B   Electronics    Model    2120MB  Power    Amplifiers    operating    in   parallel. 
The   maximum  attainable   acceleration    of    the  piston   with   this    equip- 
ment   was    about    65   g's    in    free   air    at    100   Hz.      The    frequency   and 
power    level    of    the  piston   were   controlled  by   a    General   Radio    Il6l-A 
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Coherent    Decade   Frequency    Synthesizer.       Frequencies    can   be   selected 
to   within    less    than   0.01   Hz   with    this    synthesizer.      The  piston   mo- 
tion  was    monitored  by   measuring    the   output    of   an    Endeveco  Model    2215 
Acceler ometer    mounted   within    the  piston.       Distortion    in    the  piston 
motion,    as    determined  by    the   acceler ometer ,    set    the   lower    limit    to 
which   reliable   measurements    of   the   non-linear    distortion    in    the 
pressure   waveform   could  be   made.      Measurement    of   the  pressure   wave- 
form   in   the   cavity   was    made    in    the    far    corner    from   the  piston 
source.      Refer    to   Figure   2.      The    sound  pressure    levels    and   harmonic 
components    of   the  pressure   wave    in    the   cavity   were   monitored  by   a 
1/4    inch    diameter    Bruel    and  Kjaer    type   4l3&    Condenser    Microphone. 
This    microphone   was    mounted   such    that    the  active   monitoring    element 
of   the  microphone   was    flush  with   the    interior    of    the   cavity.      The 
microphone   was    fitted  with   an  0-ring    and   seated    in    its    orifice   with 
the   same  high-vacuum  silicon   grease.      The   output    signal    of    the 
microphone   was    amplified   by   a    Bruel    and  Kjaer    type   2801   Power 
Supply.      This   power    supply    with   the   microphone    comprized   the   mi- 
crophone  system.      The   microphone   system  was    monitored   by   a    Hewlett- 
Packard  Model    400D  Vacuum   Tube   Voltmeter .      Another    vacuum    tube 
voltmeter    of    the   same   type   was   used   to   monitor    the   output    of    the 
acceler ometer .       A   Hewlett-Packard  Model    302A    Wave   Analyzer    was 
used   to  measure   harmonic    distortion    of  both    the   acceler ometer    and 
microphone   output    signals.      A  Hewlett-Packard  Model    521A    Electronic 
Counter    was    used   to  monitor    the   filtered   output    of   the   wave   ana- 
lyzer   to  assure    that   the   proper    frequency   was    being   analyzed. 
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B.   EQUIPMENT  CALIBRATION 

Absolute  calibration  of  the  microphone  system  response  was 
required  to  determine  strength  parameters.   The  microphone  system 
included  the  microphone,  its  preamplifier,  and  the  vacuum  tube 
voltmeter  connected  to  the  output  of  the  preamplifier.   The  micro- 
phone system  response  is  defined  by 

S   =  V  /P 
m     m 

where  V   is  the  RMS  voltage  from  the  preamplifier  as  indicated  on 
the  VTVM  and  P  is  the  RMS  pressure  amplitude  at  the  microphone 
face.   Calibration  of  the  microphone  system  was  accomplished  with 
the  use  of  a  Breul  and  Kjaer  Model  4220  Pistonphone  which  produces 
a  124  +,  0.2  re  volt/ubar  acoustic  signal  at  250  Hz.   The  micro- 
phone system  sensitivity  was  found  to  be 

S   =  (1.46  +  0.C3)  10~3  volt/(Nt/m2) , 

which  is  in  close  agreement  with 

S   =  1.50  +  0.03  x  10"3  Volt/(Nt/m2) 
m        — 

and 

S   =  1.52  +  0.03  x  10  J  Volt/(Nt/m  ) 
m        —  v 

as  determined  by  Beech  [2]  and  Winn  [3]  respectively.   The  micro- 
phone system  had  a  flat  frequency  response  to  within  +_  0.5  db 
from  80  Hz.  to  20  kHz.,  as  determined  from  the  microphone  response 
curve  supplied  by  the  manufacturer.   With  the  microphone  system 
calibrated,  the  wave  analyzers  were  calibrated  with  respect  to  the 
microphone  system  VTVM.   This  was  accomplished  by  removing  the 
microphone  and  its  preamplifier  from  the  VTVM,  and  substituting 
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a  voltage  signal  from  the  frequency  synthesizer  to  both  the  VTVM 

and  the  wave  analyzers  simultaneously. 

The  calibration  of  the  acceler ometer  was  accomplished  by 

direct  measurement  of  the  piston  face  displacement  at  100  Hz.  and 

at  various  power  settings  with  a  stroboscopic  light  and  a  traveling 

microscope.   The  frequency  of  the  strobe  light  was  set  very  near 

100  Hz  so  that  the  beats  of  the  piston/strobe  light  interaction 

could  be  observed.   These  beats,  which  occur  at  the  difference 

frequency   f        -  f  .      ,  allowed  the  piston  travel  to  be 
M    '  '  strobe    piston  '  ^ 

easily  measurable  with  the  traveling  microscope.   Twelve  deter- 
minations of  the  acceler ometer  sensitivity  resulted  in 

S   =  (2.36  +0.04)'  10"H  Volt/(m/sec  ). 

This  compares  with  the  value 

SA  =  (2-35  +  0.03)*  10"4  Volt/(m/sec  ) 
determined  by  Beech  L33. 

C.   MICROPHONE  OUTPUT  AND  STRENGTH  PARAMETER 

In  order  to  investigate  the  pressure  waveform  in  the  cavity, 
as  measured  at  the  microphone,  it  was  necessary  to  relate  the 
strength  parameter  to  observable  quantities.   The  strength  para- 
meter is  defined  as 

SP  =  MbQ  . 

The  peak  pressure   P,    is    related   to   the   RMS   pressure  P  by 

K  ^  1  ^  rms 


P1    =  \|T  Pr 


ms 


It  is  known  from  the  microphone  sensitivity  that 

P     =  V/S 

rms       m 
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where  V  is  the  RMS  voltage  reading  on  the  microphone  system  VTVM. 
The  strength  parameter  can  now  be  reformulated  as 

SP  =  {2    VbQ./(S  pc2) , 
1    m 

all  of  which  are  observable  or  calculable  quantities.   The  calcu- 
lation or  measurement  of  each  of  these  quantities  introduces  error 
into  the  determination  of  the  strength  parameter;  this  error  is 
estimated  to  be  about  1.0  to  1.5  percent  for  both  the  modes 
investigated. 

D.  FREQUENCY  PARAMETER 

The  frequency  parameter  is  a  quantity  which  indicates  the 
position  of  the  frequency  being  analysed  relative  to  the  fundamental 
resonance  of  the  system.   The  frequency  parameter  is  defined  by 

FP  =  2(W-Wlr)Q1A0lr 

and  can  be  shown  to  be  related  to  the  experimental  frequency 

increment  by 

FP  =  2(OJ-C0lr)/(C01  +  -601_). 

E.  SYSTEM  ALIGNMENT 

The  lower  limit  to  which  significant  measurements  could  be 
taken  was  determined  by  the  distortion  introduced  by  the  piston 
motion.   In  an  effort  to  obtain  meaningful  data  at  the  higher 
harmonics  and  for  any  harmonics  at  low  strengths,  great  care  was 
taken  to  reduce  the  distortion  introduced  into  the  cavity.   The 
distortion  of  the  piston  motion  when  free  of  the  cavity  was  below 
the  general  electronic  noise  level  of  the  acceler ometer  system  and 
thus  was  insignificant.   The  bulk  of  the  distortion  came  from 
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misalignment  of  the  piston  and  cavity  orifice.   The  procedure  used 
to  obtain  maximum  alignment  of  the  piston  and  thus  minimize  the 
distortion  was  as  follows:  The  output  signal  of  the  acceler ometer 
system,  which  measured  the  piston  motion,  was  displayed  on  the 
dual  trace  oscilloscope  and  this  waveform  was  compared  to  a  pure 
sine  wave.   The  piston  and  exciter  apparatus,  which  was  mounted  in 
an  adjustable  framework,  was  positioned  relative  to  the  cavity 
orifice  until  the  acceler ometer  waveform  was  as  clean  and  undis- 
torted  as  possible.   In  order  to  minimize  any  effects  which  may 
have  arisen  from  vibration  coupling  of  the  cavity  and  exciter  or 
cavity  and  table,  the  cavity  was  mechanically  isolated  from  both 
the  exciter  and  the  table.   This  was  accomplished  by  mounting  both 
the  cavity  and  the  piston-exciter  apparatus  on  platforms  con- 
structed of  alternate  slabs  of  wood  and  rubber.   This  restricted 
the  interaction  of  the  cavity  with  the  piston  to  the  piston  ori- 
fice only.   With  proper  alignment,  it  was  possible  to  reproduce 
conditions  where  the  distortion  in  the  input  waveform  to  the 
cavity  was  less  than  0.2%  second  harmonic  and  0.3%  for  third  har- 
monic, as  measured  by  the  H-P  302A  wave  analyzer.   Distortion 
values  for  the  higher  harmonics  were  below  the  general  ambient 
noise  of  the  total  apparatus,  and  thus  were  ignored. 
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IV.   PROCEDURES 

A.   PRERUN  PROCEDURES 

Before  any  meaningful  data  could  be  taken,  it  was  important  to 
determine  the  amount  of  vibration  transferred  through  the  cavity 
walls  and  picked  up  by  the  microphone.   To  make  this  determination, 
the  top  of  the  cavity  was  removed,  thus  relieving  the  pressure  seal, 
and  a  wood  plate  added  between  the  piston  and  microphone  to  isolate 
the  microphone  from  the  pressure  waves  induced  in  the  air  by  the 
piston.   As  described  in  the  System  Alignment  section,  the  cavity 
and  the  exciter -piston  assemblies  were  isolated  by  a  floating 
platform  technique.   The  vibration  picked  up  by  the  microphone  was 
measured  at  frequencies  of  interest  and  produced  microphone  system 
responses  of  8*10  '  to  2*10    volts.   In  all  cases  these  values 
amounted  to  distortion  well  below  0.3%,  which  is  the  maximum 
amount  of  piston  motion  distortion. 

It  was  noted  by  Beech  [3]  that  temperature  changes  produced 
significant  changes  in  the  resonant  frequency  of  the  tube  he  in- 
vestigated.  The  temperature  effect  on  the  cavity  investigated 
here  was  very  pronounced.   This  presented  many  problems  as  the 
frequency  of  fundamental  resonance  was  the  reference  for  all  de- 
terminations of  frequency  parameter.   Therefore  it  was  necessary 
to  know  the  value  of  the  fundamental  resonance  at  all  times.   To 
minimize  this  shift  in  the  resonance  frequency,  the  entire  ap- 
paratus was  allowed  to  warm  up  for  a  period  of  one  to  four  hours, 
depending  on  the  strength  parameter  to  be  investigated,  prior  to 
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taking  data.   This  procedure  did  not  totally  relieve  the  effect  of 
the  frequency  shift  so  another  calibration  method  was  employed. 
Many  determinations  of  the  fundamental  frequency  were  made  during 
the  taking  of  data,  thus  producing  a  calibration  curve,  see 
Figure  3>  by  which  the  frequency  parameter  values  could  be 
corrected. 

Immediately  prior  to  the  taking  of  each  data  set,  the  align- 
ment of  the  piston  and  cavity  was  checked  as  outlined  in  the  system 
alignment  section  of  Experimental  Considerations.   From  this  point 
onward,  the  procedure  for  data  collection  was  dictated  by  the 
quantities  being  measured.   There  was  one  experimental  method  to 
measure  the  Q's  and  E's  of  the  system  in  the  linear  region,  and 
another  to  measure  the  harmonic  content  of  the  non-linear  finite 
amplitude  waves. 

B.   MEASUREMENT  OF  Q(N)  AND  E(N)  IN  THE  LINEAR  REGION 

Data  collection  in  this  region  required  very  little  power, 
but  required  very  fine  control  of  the  power  setting.   This  control 
was  obtained  by  inserting  a  voltage  dividing  circuit  between  the 
frequency  synthesizer  and  the  power  amplifiers.   This  voltage  di- 
vider reduced  the  available  power  by  a  factor  of  ten  and  increased 
the  sensitivity  of  the  power  control  by  the  same  factor  of  ten. 
With  this  arrangement,  a  satisfactory  fine  tuning  power  control 
was  obtained  at  low  strengths. 

With  the  power  set  at  some  nominal  low  value,  the  frequency 
region  about  the  resonant  or  harmonic  pressure  peak  of  interest 
was  scanned  until  a  maximum  in  pressure  was  indicated  by  the 
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microphone  system  response.   The  input  power,  determined  in  the 
linear,  low  strength,  region,  and  kept  constant  during  the  fol- 
lowing measurements.   The  frequency  of  the  driving  piston  was 
changed  above  and  below  the  frequency  of  maximum  pressure,  until 
the  pressure  was  -3  dB  of  the  maximum  pressure  value.   This  deter- 
mination of  the  -3  dB  points  was  made  with  the  302A  wave  analyzer , 
which  measured  the  pressure  amplitude  of  the  cavity  by  way  of  the 
microphone  system. 

It  was  found  that  the  first  seven  overtones  of  the  nOO  mode 
and  the  first  five  overtones  for  the  OnO  mode  were  undisturbed 
by  overlapping  Q  curves.   The  value  of  Q  for  each  response  curve 
was  determined  from  the  frequency  difference  of  the  -3  dB  points, 
divided  by  the  center  point  frequency  between  the  two  -3  dB 
points.   This  center  point  was  assumed  to  be  a  good  estimate  of 
the  peak  pressure  frequency,  CO   ,  for  the  response  curve  in 
question.   It  was  with  these  values  for  the  system  resonance  fre- 
quencies that  the  E(n)'s  were  determined. 

C.   MEASUREMENT  OF  THE  HARMONIC  CONTENT  OF  FINITE  AMPLITUDE 

STANDING  WAVES 

In  order  to  obtain  finite  amplitude  waves  in  the  cavity,  it 
was  necessary  to  go  to  much  higher  strengths  than  were  required 
for  the  measurement  of  the  Q's  and  E's.   With  the  driving  fre- 
quency set  to  some  value  near  the  fundamental,  and  the  peak  pres- 
sure component  of  the  fundamental  held  constant,  the  strength  of 
each  harmonic  was  measured  by  one  of  the  wave  analyzers.   The 
other  wave  analyzer  monitored  the  fundamental  pressure  P  ,  which 
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was  kept  constant  by  varying  the  power  setting  of  the  synthesizer 
as  necessary.   Values  of  the  frequency  parameter  were  chosen  such 
that  the  distortion  curves,  Figures  4  through  9,  could  be  deter- 
mined in  the  areas  of  maximum  distortion.   The  order  in  which 
measurements  were  taken  for  each  distortion  curve  was  randomized 
to  avoid  systematic  error.   Strength  parameters  of  about  0.5,  0.75, 
and  1.0  were  investigated  for  each  mode.   Strengths  below  0.5  did 
not  produce  high  enough  distortion  to  be  measured  with  the  equip- 
ment available.   Strengths  above  1.0  were  not  accessable  due  to 
limitations  in  power . 
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V.   RESULTS 

Tables  I  and  II  are  the  data  for  the  linear,  infinitesimal 
amplitude  region.   The  information  on  the  empirical  losses  is  con- 
tained in  the  Q's  and  E's.   These  are  the  values  used  in  the  com- 
puter program  to  predict  the  harmonic  distortion  curves  in  Figures 
4  through  9.   The  theoretical  curves  in  these  figures  are  plotted 
along  with  the  experimentally  measured  values  of  harmonic  dis- 
tortion.  The  theoretical  predicted  values  were  generated  for  fre- 
quency parameter  intervals  of  0.2  and  were  plotted  as  solid  curves. 
The  experimentally  measured  values  were  plotted  as  dots  where  they 
occured.   Data  were  taken,  and  theoretical  predictions  made,  for 
strength  parameters  of  0.513,  0.782,  and  1.00  for  the  nOO  mode  and 
0.5Q7,  0.780,  and  1.00  for  the  OnO  mode. 

The  point  of  intersection  of  the  error  bars  on  the  P   /P      theo- 
retical curve  in  Figure  4  corresponds  to  the  value  of  EP  =  2Q  E(2) 

60 
2r 
where  40  =  — - —  .   It  is  important  to  note  that  the  n  =  2  distortion 

curves  peak  when  the  system  is  driven  at  this  frequency.   That  is, 

when  the  driving  frequency  CO  is  equal  to  the  classical  fundamental 

of.CO^  ,  there  is  maximum  content  of  P„.   The  point  where  this  oc- 
2r  2        t 

curs  for  each  P   /P      curve  is  indicated  by  arrows  in  Figures  5 

through  9  and  by  crossed  error  bars  in  Figure  4.   The  error  in  the 
determination  of  these  points  comes  mainly  from  the  error  in  E(2) 

that  appears  in  the  above  equation  for  FP .   The  error  in  Q   is  of 

much  less  significance  as  is  indicated  in  Table  III.   The  error  in 

the  P/P  curves  in  Figures  5  through  9  are  approximately  the  same 
as  those  indicated  in  Figure  4. 
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Table   1 


nOO 

Mode 

n 

f_(Hz) 

f+(Hz) 

Af (Hz) 

f   (Hz) 

n          ' 

Q 

E  x    103 

1 

677.80 

679.70 

1.90 

678.75 

357.24 

0 

2 

1365.74 

1368.49 

2.75 

1367.12 

497.13 

7.0 

3 

2059.17 

2062.84 

3.67 

2061.01 

561.58 

12.0 

4 

2787.37 

2792.73 

5.36 

2790.05 

520.53 

27.5 

5 

3423.00 

3427.80 

4.80 

342  5.40 

713.63 

9.1 

6 

4126.41 

4131.91 

5.50 

4129.16 

750.76 

13.6 

7 

4803.92 

4809.93 

6.01 

4806.93 

799.82 

H.4 

1 

678.02 

679.95 

1.93 

678.99 

351.81 

0 

Table    II 


OnO 

Mode 

n 

fjHz) 

f+(Hz) 

Af(Hz) 

f    (Hz) 

Q 

E   x    103 

1 

813.99 

816.08 

2.09 

815.04 

389.97 

0 

2 

1633.16 

1636.39 

3.23 

1634.78 

506.12 

2.8 

3 

2415.11 

2420.57 

5.46 

2417.84 

442.83 

-11.3 

4 

3143.74 

3150.10 

6.36 

3146.92 

494.80 

-34.9 

5 

3911.05 

3920.35 

9.30 

3915.70 

423.62 

-39.4 

1 

814.23 

816.34 

2.11 

815.29 

386.39 

0 
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Frequency   Parameter    (FP) 
nOO   Mode  SP    =   0.513 

Figure   4 
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Figure    5 
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Frequency   Parameter     (FP) 
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Figure   8 
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0.002 


Frequency    Parameter     (FP) 
OnO   Mode  SP    =    1.00 

Figure   9 
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Table  III  contains  the  results  of  five  separate  determinations 
of  the  Q's  and  E's  taken  in  the  infinitesimal  amplitude  region. 
These  measurements  were  made  to  provide  an  estimate  of  the  error 
in  the  theoretical  curves.   In  all  cases  the  values  used  to  gene- 
rate the  curves  were  very  close  to  the  average  values  determined 
later  and  the  regeneration  of  the  curves  was  therefore  unwarranted, 
Table  III  contains  both  sets  of  values  for  direct  comparison. 
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nOO   Mode 


Table   III 


Values    Used   to 
Average   Values  Determine   Curves 


Q(n) 


1  358.2  +  0.4 

2  498.5  +  4.7 

3  559.8  +  0.8 

OnO  Mode 

1  359.8   +  5.3  

2  515.7  +  5.1  2.6  +  0.1 

3  459.2  +  7.5       -11.3  +  o.l 


E(n)*103 

Q(n) 

E(n).103 

357.2 

6.8   +0.1 

497.1 

7.0 

11.6    +   0.2 

561.6 

12.0 

390.0 

509.1 

2.8 

442.8 

-11.3 

33 


VI.   CONCLUSIONS 

It  appears  that  the  present  theoretical  model  successfully 
predicts  the  major  features  of  the  harmonic  content  for  finite- 
amplitude  standing  waves  in  a  rigid-walled  cavity.   In  all  cases 
investigated,  the  shape  and  general  magnitude  were  predicted;  see 
Figures  4  through  9. 

It  is  interesting  to  note  the  position  of  each  theoretical 
curve  relative  to  its  corresponding  experimental  curve.   For  the 
nOO  mode,  the  theoretically  predicted  values  are  consistently 
above  the  measured  values.   The  situation  for  the  OnO  mode  is  the 
same  for  P   /P      but  opposite  for  P^/P, .   A  possible  source  of  this 
difference  is  the  large  difference  in  the  geometry  encountered  by 
each  mode.   The  nOO  mode  is  a  standing  wave  whose  wave  vector  is 
perpendicular  to  the  wall  with  the  microphone  orifice,  which  must 
be  considered  a  very  small  perturbation  due  to  its  small  size.   On 
the  other  hand,  the  OnO  mode  is  a  standing  wave  with  its  wave 
vector  perpendicular  to  the  piston  orifice,  which  is  a  large  per- 
turbation in  the  geometry.   In  addition,  the  fact  that  the  wave 
generating  surface,  the  piston  face,  is  normal  to  the  nOO  wave 
vector  and  perpendicular  to  the  OnO  vector  could  cause  a  con- 
sistant  difference  in  the  behavior  of  the  modes.   It  is  certainly 
unclear  to  this  investigator  just  what  the  specifics  are  that 
explain  the  trends  in  the  theory/experiment  relationship.   It  is 
clear  however,  that  the  theory  of  Coppens  and  Sanders  as  modified 
to  include  empirical  losses,  does  predict  the  essence  of  the 
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behavior  of  finite-amplitude  standing  waves  in  the  rigid-walled 
cavity  investigated. 
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APPENDIX  A 

In  the  investigation  by  Winn  [2],  it  was  found  that  the  tubes 
analyzed  did  not  behave  in  a  manner  consistent  with  Rayleigh- 
Kirchoff  wall  losses.   In  an  effort  to  determine  if  the  interior 
finish  of  the  tubes  was  a  critical  factor  that  determined  the  loss 
mechanisms,  the  six  foot  tube  was  reanalyzed  after  refinishing  the 
interior  to  a  highly  polished  state.   These  results  were  compared 
to  the  results  for  the  tube  in  a  smooth  but  not  polished  condition 
as  determined  by  Winn. 

It  was  found  that  the  absorption  coefficients,  defined  as 

CO   -  CO 

a     =     n+   n" 


n      2C 
o 

were  not  significantly  changed  from  the  values  found  by  Winn,  and 
both  sets  of  values  obeyed  the  Rayleigh-Kir choff  predictions 
indicated  by  the  solid  line  in  Figure  10.   The  behavior  of  the 
overtones  (for  the  polished  tube)  was  found  to  be  significantly 
different  (from  the  overtone  behavior  in  the  unpolished  tube), 
and  was  significantly  different  from  the  Rayleigh-Kirchof f  pre- 
dictions as  indicated  by  a  solid  line  in  Figure  11. 
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